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Measurements of time dependent CP asymmetry in B ^ VV decays with 
BELLE 

Ryosuke Itoh^ 

Representing the Belle Collaboration 

'^Institute of Particle and Nuclear Studies, High Energy Accelerator Research Organization (KEK), 
1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 

A study of CP violation in _B — > J/tpK* (Kg-K^) decays by time dependent angular analysis is discussed. Status 
of time independent analyses for other B —> VV decays is also reported. The data used for the analyses are taken 
with the Belle detector at KEK. 



1. Introduction 

There are three helicity states in B ^ VV de- 
cays. Although one of the states is a pure CP even 
state, CP even and odd states are mixed in other 
states and the observed CP asymmetry is diluted. 
By studying decay angles of final state particles, 
it is possible to project out each CP state in a 
statistical way. A theoretical expression of differ- 
ential time-angular decay rate of B — > VV decays 
consists of components corresponding to three he- 
licity states and also to the interference between 
them. Each term is expressed as a product of an 
angular term and an amplitude term. The am- 
plitude term contains CP violating phase(s) as 
a function of At, the decay time difference be- 
tween two B mesons from an T{AS) decay. A fit 
of the expression for the measured decay angles 
and At gives the determination of the CP violat- 
ing phase(s). The interference terms are also rich 
sources of various interesting physics such as the 
cos 2(^1 measurement in B° J/ipK*°{Kg7r") 
decays, and the simultaneous determination of r 
and sin(2(/)i -I- (p^) in B° D*~ p'^ decays where 
r is the ratio of doubly-Cabibbo-suppressed mode 
to Cabibbo favored mode contributing to the de- 
cay. 

In this talk, a study of time dependent an- 
gular analysis for B J/^K*{KgTr'^) is dis- 
cussed. Also reported is status of angular anal- 
yses for B" D*+D*- and B° D*-p+. 
The data used in the analyses are taken with 



the Belle detector in the KEKB accelerator^. 
The KEKB accelerator consists of two separate 
rings for electrons and positrons. The energies 
are set at 3.5 GeV for positrons and 8.0 GeV 
for electrons, respectively. These two beams 
are collided at the interaction point to produce 
T(4S') in motion for the measurement of At. 
The peak luminosity of the machine has reached 
7.4 X 10^^cm~^ sec~^ which is the current world 
record. The total integrated luminosity up to now 
is 89.1/6"^. The Belle detector is placed at the 
interaction point of the KEKB accelerator. It is 
a general-purpose detector with a wide angle cov- 
erage featuring a particle identification using the 
Aerogel Cherenkov Counter assembly. The data 
sample used in the analyses corresponds to an in- 
tegrated luminosity of 78 fb~^ taken on the T(45') 
resonance. 

2. Time dependent angular analysis for 

Candidate mesons are reconstructed by se- 
lecting events with a J/^ identified from a pair of 
oppositely-charged leptons and a K* from a pair 
of Kg and 7r° candidates. The beam-constrained 
mass (Afbc), which is the invariant mass of a re- 
constructed J/4> and K* calculated taking the 
energy to be the beam energy, is required to be 
in the range 5.27 — 5.29 GeV/c^ and the energy 
difference between B candidate and the beam en- 
ergy (AE) to satisfy -50 MeV < A£; < 30 MeV. 
To eliminate slow 7r° backgrounds, the angle of 



2 



the kaon with respect to the K* direction in 
the K* rest frame, cos 9k* , is required to satisfy 
cosOk* < 0.8. When an event contains more than 
one candidate passing the above requirements, 
the combination for which AE is closest to zero 
is selected. The number of events remaining after 
the selection is 104. 

The time dependent angular distribution of the 
decay is described asQ: 

d''r{0trAtr,9K',At) 



d cos 9trd(f>trd COS 9k' dAt 

p-|At|/r^o 6 



2t 
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J29^{^trAtr,9K')a^{At) (1) 
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where gi are the angular terms and are ampli- 
tude terms, tb is the lifetime of a meson. The 
gi are expressed as: 



gi = 2cos^ 6'i<--(l - sin^ cos^ (2) 
sin'^ 6k* (1 — sin^ 6tr sin^ (j)tr) (3) 

(4) 
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g^ ~ sm^ Ok* sin^dtr 

-1 . . 2 

54 = — p sin 29k* sin 9tr sin 2(/)tr 
v2 

(75 sin'^ 9k* sin 20^^ sin 

ge = —z= sin 26'if » sin 29tr cos (/)tr 
v2 



(5) 
(6) 
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where 0tr, (t>tr and are the decay angles de- 
fined in the transversity basis The are 

ai = |AoP(l + 77sin2(/)i sin AmAt) (8) 

a2 = |A|||^(l-|-rysin20isinAmAt) (9) 

flg = |A_Lp(l - 77sin20i sin AmAi) (10) 
04 = i?e(A|Ao)(l + 77sin20i sin AmAi) (11) 
a5 = ry/m(A|^^) cos AmAi — 

r7i?e(^|^j_)cos2(/)i sinAmAt (12) 
ae = 'r]Im{AQA±) cos AmAt — 

r]Re{A'^A±_) cos2(j)i sin AmAt (13) 

where Am is the B'^ — mixing parameter, t] 
is -1-1 for B° while -1 for W. Aq, Ay and A± 
are the complex decay amplitudes for three he- 
licity states in the transversity basis. Two CP 
violation parameters, sui2(f>i and cos20i, appear 
in the formula. The values of these parameters 



are determined by fitting the formula to the mea- 
sured angles and At, taking into account the de- 
tection efficiency and background. The fit is done 
using an unbinned maximum likelihood method. 
The probability density function for an event is 
defined as 

PDF = Ug{Ml,rMGtr,<ptr,9K*) 
d'*r(cOS 6tr, 4>tr, COS 9k* , At) 



dcOS Otrdijitrd COS 9k* dAt 
-1 At I /t o 



+ 



X {J2 fcf{Mbc)ADF^f{9tr, <l>tr, 9 k*) 

i 

+ fnr{Mbc)ADFnr{9tr,(t>tr,9K*)} 

+ S{At)f,,^u{Mt^)ADF,,^u{9tr, (ptr, 9k*). (14) 

fsig, fcf, fnr and fcombi are fractions of 
the signal, cross feed, non-resonant production, 
and combinatorial background as a function of 
Mfcc, respectively, while AD Fcf, AD Fnr and 
ADFcombi are corresponding angular shape func- 
tions, e is a three dimensional detection efficiency 
function for the signal. The determinations of 
these functions and the decay amplitudes are de- 
scribed in ref.|^. Phases in the amplitudes are 
chosen so as to conserve s-quark helicity||^. The 
Tb and Am are fixed at the PDG values in the fit. 
The flavor tagging procedure [|| gives the flavor of 
tag-side B meson q and the probability of wrong 
tag w. To account for the effect of the wrong tag- 
ging in the data, ij is replaced with —q{l — 2w) in 
the PDF. 

Each term in the PDF is then convolved with 
appropriate resolution functions |^ separately for 
the signal, the backgrounds with a B'^ lifetime 
and the combinatorial background with a 5{At) 
function shape. The parameters in the resolution 
functions are calculated event by event. From 
the fit to the data, we obtain sin2</)i = 0.13 ± 
0.51 ±0.06, and cos2(/)i = 1.40 ± 1.28 ±0.19. The 
systematic errors include uncertainties in the res- 
olution parameters, wrong tagging fractions, de- 
cay amplitudes, background fractions and shapes, 
and others. The distributions of At measured for 
the samples tagged as g = +1 and —1 are shown 
in Figure |l| together with the predictions by the 
PDF with obtained CP parameter values. 
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Figure 1. At distributions for samples tagged as 
q — +1 and —1. The lines show the predictions by 
the PDF with determined CP parameter values. 



Figure 2. Distribution of cos 0tr for B' 
D*+D*~ candidates. 



3. Time independent angular analyses 

The time dependent angular analysis is appli- 
cable for other B VV decays such as i?" 
D*+D*-, D-*p+ and for the determi- 

nation of CP violating phases including those 
other than 0i . Time dependent analyses for these 
modes are still in preparation and the status of 
time independent (integrated) analyses is reported 
here. The measurement of the decay amplitudes 
by the time independent angular analysis is es- 
sential for the determinations of the CP violating 
phases. 

3.1. B" D*+D*- 

The reconstruction of a D* is performed us- 
ing a Dtt decay with a subsequent D decay into 
Kit, Kinr, JCtttttt, KK and KKt:. A D meson is 
identified when the invariant mass of decay par- 
ticles is within 3 to 6 a from the nominal D 
mass depending on the decay mode. A D* is 
reconstructed by calculating the mass difference 
M^Dtt) — M{D) and those satisfy the difference 
within S.OMeFpO) or 2.2^MeV(D+) are iden- 
tified as D* . The candidate B mesons are se- 
lected by applying M^c > Mga — 3cr and |Ai?| < 
AOMeV. 

The distribution of the transversity angle |^ is 
measured in the selected B^ candidates as shown 



in Fig. g. The distribution can be expressed as 
a function of the fraction of CP-odd component 

{Rt) as 

— = 7(1 - -Rt) sin'' Otr + \Rt cos" 9tr. (15) 

The predictions given by this formula for cases 
with Rt = and Rt = 1 are also shown in the 
figure. As seen, the distribution of the data is 
close to the prediction with Rt = 0. 

3.2. B° D*-p+ 

The reconstruction of a D*^ is performed in a 
similar manner as described above. A p+ is re- 
constructed from a pair of 7r+ and tt*^ candidates. 
The Mbc of a pair of reconstructed D* and p is 
required to be in the range 5.27 - 5.29 GeV/c^ to 
and A-B to satisfy -100 < \l^E\ < 50MeV. 

The time independent angular distribution is 
studied using the helicity basis angles The de- 
cay rates as a function of each of three projected 
angles (two helicity angles 9i where i = 1,2, and 
an angle between two decay planes x) can be ex- 
pressed as 

= cos e. + -{\A^\ +\A,\ )sm e., 

(16) 

.4 |2_^ 8 2. . 2 
^^-(gi^ol +glAJ )sm X 



^Thnangular analysis for _B+ — > p^p" is covered in other 
talkjig. 



+ (^i^oi' + l\A^^\')cos'x+llm{AlA^)sm2x 

(17) 
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Figure 3. Distributions of angles used in the an- 
gular analysis for B'^ D*~ candidates. The 
distributions are corrected for the acceptance and 
background contaminations. 



The decay angles in the transversity basis {9tr, 
4>tr and Ok*) and the decay time difference of 
two B mesons (Af) are measured for each event 
with the determination of the flavor of tag-side 
B meson. The values are fitted to the theoretical 
distribution using an unbinned maximum likeli- 
hood method considering the effects of detector 
acceptance, background contamination, At reso- 
lution and wrong flavor tagging. From the fit, 
two CP violation parameters are determined to 
be sin 201 = -h0.13 ± 0.51 ± 0.06, and cos20i = 
+l.m ± 1.28 ± 0.19. If we take the s-quark he- 
licity conservation choice of the amplitude phases 
and the sin20i value measured using other de- 
cay modes (0.72 ± 0.07 ± 0.04) @, the obtained 
sign of cos 201 suggests the 201 to be in the range 
0° < 201 < 90° and the other choice of 20i > 90° 
is not preferred although the current statistics is 
not enough to conclude this. 

A time independent angular analysis for other 
decay modes, B° D*+D*-, B° D*-p+, and 
B~^ —> p'^ , are now being performed. 



where Aq, A\\ and A±_ are the decay amplitudes 
in the transversity basis. In addition, we also used 
X-differential decay rates projected in 4 quadrants 
of a 01 — 62 plane (angular regions from to 7r/2 
and 7r/2 to tt for each of them) to observe the 
effects in interference terms. The difference in 
the decay rates between two diagonally-summed 
quadrants (A) can be written as 



dx 
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Im{A*)AA_) sin x H 7=-Re(^o^|| ) cos x- 

9\/2 

(18) 



Fig. m shows the distributions of cos 6i {61 and 
6*2 distributions are added), x^ and A(x). The 
measurement of the decay amplitudes is now in 
progress by fitting these distributions with the 
theoretical formula given above considering ef- 
fects of the detector acceptance and background 
contaminations. 

4. Summary 

A full time dependent angular analysis is per- 
formed for B° J/iI:K*^{KItt^) decays col- 
lected with the Belle detector at KEK B-factory. 
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